Light Scalar Meson er(600) in QCD Sum Rule with Continuum 
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The light scalar meson <t(600) is known to appear at low excitation energy with very large width 
on top of continuum states. We investigate it in the QCD sum rule as an example of resonance 
structures appearing above the corresponding thresholds. We use all the possible local tetraquark 
currents by taking linear combinations of five independent local ones. We ought to consider the 
n-TT continuum contribution in the phenomenological side of the QCD sum rule in order to obtain a 
good sum rule signal. We study the stability of the extracted mass against the Borel mass and the 
threshold value and find the cr(600) mass at 530 MeV ± 40 MeV. In addition we find the extracted 
mass has an increasing tendency with the Borel mass, which is interpreted as caused by the width 
of the resonance. 
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I. INTRODUCTION 

The light scalar mesons, cr(600), k(800), /o(980) and ao(980), have been intensively discussed for many years PHI]. 
However, their nature is still not fully understood 041] ■ They have the same quantum numbers J PC = ++ as the 
vacuum, and hence the structure of these states is a very important subject in order to understand non-perturbative 
properties of the QCD vacuum such as spontaneous chiral symmetry breaking. They compose of the flavor SU(3) nonet 
with the mass below 1 GcV, and have a mass ordering which is difficult to be explained by using a qq configuration 
in the conventional quark model [9hfl3j . Therefore, several different pictures have been proposed, such as tetraquark 
states and meson-meson bound states, etc. Here we note that hadrons with complex structures such as tetraquarks 
may exist in the continuum above the threshold energy of two hadrons with simple quark structure. 

The tetraquark structure of the scalar mesons was proposed long time ago by Jaffe with an assumption of strong 
diquark correlations [HI, HH • It can naturally explain their mass ordering and decay properties ■ Yet the basic 

assumption of diquark correlation is not fully established. In this letter, we study er(600) as a tetraquark state in the 
QCD sum rule approach as an example of resonances in the continuum states above the 7r-7r threshold. In the QCD 
sum rule, we calculate matrix elements from the QCD (OPE) and relate them to obscrvables by using dispersion 
relations. Under suitable assumptions, the QCD sum rule has proven to be a very powerful and successful non- 
perturbative method in the past decades (l9l . |20| . Recently, this method has been applied to the study of tetraquarks 
by many authors plM^ . 

In our previous paper J24[, we have found that the QCD sum rule analysis with tetraquark currents implies the 
masses of scalar mesons in the region of 600 - 1000 MeV with the ordering m a < m K < mf oMo , while the conventional 
qq current is considerably heavier (larger than 1 GcV). To get this result, first we find there arc five independent local 
tetraquark currents, and then we use one of these currents or linear combinations of two currents to perform the QCD 
sum rule analysis. But these interpolating currents do not describe the full space of tetraquark currents. In order 
to complete our previous study, we use more general currents by taking linear combinations of all these currents. It 
describes the full space of local tetraquark currents which can couple to a (600). Since ct(600) meson is closely related 
to the 7T-7T continuum and it has a wide decay width, we also consider the contribution of the 7r-7r continuum as well 
as the effect of the finite decay width. 

This paper is organized as follows. In Sec. UH we establish five independent local tetraquark currents, and perform 
a QCD sum rule analysis by using linear combinations of five single currents. In Sec. IIII1 we perform a numerical 
analysis, and we also study the contribution of 7T-7T continuum. In Sec. IIV1 we consider the effect of the finite decay 
width. Sec. [V] is devoted to summary. 
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II. QCD SUM RULE 



The local tetraquark currents for ct(600) have been worked out in Ref J24[. There are two types of currents: 
diquark-antidiquark currents (qq)(qq) and meson- meson currents (qq)(qq). These two constructions can be proved to 
be equivalent, and they can both describe the full space of local tetraquark currents [24| . Therefore we shall just use 
the first ones. Since we use their linear combinations to perform the QCD sum rule analysis, we can not distinguish 
whether it is a diquark-antidiquark state or a meson-meson bound state. However, we find that tetraquark currents 
with a single term do not lead to a reliable QCD sum rule result which means that cr(600) probably has a complicated 
structure. The five independent local currents are given by: 

S3 = (ulC-f 5 d b )(u a -f 5 Cd% - u b j 5 CdJ) , 

Vi = KC W d b )(u a ^-y 5 CdJ -u b7 »-y 5 C%) , 

T° = {ulCa^d b ){u a a^Cd^ + u h a» v CdF a ) , (1) 

Al = {u T a C ltl d b ){u a ^CdJ +u b ^CdF a ), 

p% = (ulcd b )(u a cdj - u b cd^) . 

The summation is taken over repeated indices v, ■ ■ ■ for Dirac, and a, b, ■ ■ ■ for color indices). The currents 5*, 
V, T, A and P are constructed by scalar, vector, tensor, axial-vector, pseudoscalar diquark and antidiquark fields, 
respectively. The subscripts 3 and 6 show that the diquarks (antidiquarks) are combined into the color representations, 
3 C and 6 C (3 C and 6 C ), respectively. 

These five diquark-antidiquark currents (qq)(qq) are independent. In this work we use general currents by taking 
linear combinations of these five currents: 

V = tie i6l Si +t 2 e i62 V 3 a + t 3 e ies Tg + Ue i0l A% +t 5 e ie5 Pg , (2) 

where ti and 8i are ten mixing parameters, whose linear combination describes the full space of local currents which 
can couple to a (600). We can not determine them in advance and therefore we choose them randomly for the study 
of the QCD sum rule. 

By using the current in Eq. ([2]), we calculate the OPE up to dimension eight. To simplify our calculation, we 
neglect several condensates, such as (g 3 G 3 ), etc., and we do not consider the a s correction, such as g 2 (qq) 2 , etc. The 
obtained OPE are shown in the following. We find that most of the crossing terms are not important such as P13, 
and even more some of them disappear: pi$ = 0, etc. For the most cases, we find that the OPE terms of Dim=6 and 
Dim=8 give major contributions in the OPE series in our region of interest. This is because the condensates (qq) 2 
(D=6) and (qq) (gqaGq) (D=8) are much larger than others. 

Since the OPE series should be convergent to give a reliable QCD sum rule, we also calculate the OPE of Dim=10 
and Dim=12. However, we find that these terms are not important. Using the parameter set (2) and the the values 
of the condensates of the next section as an example, we show the convergence of the two-point correlation function 
n(Ms, ,s ) = Jq S ° p(s)e~ s / M Bds in Fig. Q] as functions of M%. The threshold value is taken to be s = 1 GcV 2 , and 
we show its behavior up to certain dimensions. We find that the OPE up to Dim=0 and Dim=2 are very small; the 
OPE of Dim=4 gives a minor contribution; the OPE of Dim=6 and Dim=8 are both important; the OPE of Dim=10 
and Dim=12 are both small, and so we shall neglect them in the following analysis. 




Borel Mass z [GeV 1 



FIG. 1: The convergence of the two-point correlation function II(Ms,so). The threshold value is taken to be so = 1 GeV 2 , 
and we show its behavior up to certain dimensions, as functions of M%. The solid line is for H(Mb,so) up to Dim=8. The 
short-dashed line around it is for n(Ms, so) up to Dim=10, and the long-dashed line around it is for II(Mb, so) up to Dim=12. 



p{s) = t\p n {s) + t\p 22 (s) + tlp^(s) + t\p ii (s) +tlp 55 (s) 

+2* x t 2 cos (6>i - e 2 ) Pl2 {s) + 2*1*3 cos (6»i - 6 3 )p 13 (s) + 2* x * 4 cos (6 1 - 6 A )p XA {s) 
+2*2*3 cos (6*2 - 6 , 3)p23(s) + 2*2*4 cos (6 2 - 9a)P2a{s) + 2* 2 *5 cos (6» 2 - 5 )p 25 (s) 

+2*3*4 COS (6*3 - 4 )P34(s) + 2* 3 t 5 COS (6> 3 - 6» 5 )j035(s) , 



where 



.s 4 ??^ 2 m»mj _ m rf 2 3 ,(g^GG)_ _ m u (qq) _ m d (gg) 2 

PlUSj ~~ 61440tt 6 1 1536tt 6 1536tt 6 15367r ejS [ 6144tt 6 192tt 4 192tt 4 jS 

mj(g 2 GG) m u m d (g 2 GG) _ m 2 d (g 2 GG) _ m u (gqaGq) _ m d (ggaGg) (gg) 2 
1 1024tt 6 1024tt 6 1024tt 6 64tt 4 64tt 4 12tt 2 jS 

7ml(qq) 2 m u m d (qq) 2 _ 7m 2 , (qq) 2 _ m u (g 2 GG) (qq) _ m d (g 2 GG) (qq) (qq) (gqaGq) 
48tt 2 " 4tt 2 48tt 2 768tt 4 " 768tt 4 12tt 2 

s 4 m. u 2 m. u m d m d 2 3 Ag 2 GG) m u (gg) m d (qq) 2 

P22lsj - 153607r 6 + ^ 3847r 6 76g7r 6 3 847r 6jS 1 3072tt 6 24tt 4 24tt 4 JS 

m 2 u (g 2 GG) m u m d (g 2 GG) _ m 2 (g 2 GG) m u (gqaGq) m d (gqaGq) _ (qq) 2 
[ 512tt 6 512tt 6 512tt 6 32tt 4 32tt 4 6tt 2 jS 

11m 2 (gg) 2 2m u m d (qq} 2 11m 2 (gg) 2 _ m u (g 2 GG)(qq) _ m d (g 2 GG)(qq) _ (qq) (gqaGq) 
12tt 2 tt 2 12tt 2 384tt 4 384tt 4 6tt 2 

s 4 m tl 2 _ m d 2 3 ll(ff 2 GG) m u (qq) m d (qq) 2 

P33lSj ~ 1280^ 6 +l 32tt 6 3 27r 6 JS +[ 768tt6 + 4tt 4 + 4tt 4 >S 
llml(g 2 GG) _ llm 2 (g 2 GG) 
{ 128tt 6 128tt 6 ,S 

5m 2 (gg) 2 20m M m d (gg) 2 5m 2 (gg) 2 \\m u (g 2 GG) (qq) llm d (g 2 GG)(qq) 

7T 2 TT 2 TT 2 967T 4 96tT 4 

s 4 . m„ 2 m M m d _ m d 2 3 5(g 2 GG) 2 

P44lSj - 76g07r 6+l 1927r6 + gg^e 1927r 6> S + 3 07 2^6 S 

bm 2 u (g 2 GG) 5m u m d (g 2 GG) _ bm 2 d (g 2 GG) _ m u (gqaGq) _ m d (gqaGq) (qq) 2 
1 512tt 6 512tt 6 512tt 6 16tt 4 16tt 4 3tt 2 )S 

m 2 u (gq) 2 8m u m d (qq) 2 _ m 2 (gg) 2 m u (g 2 GG)(qq) m d (g 2 GG)(qq) (qq) (gqaGq) 

6tt 2 3tt 2 6tt 2 128tt 4 128tt 4 3tt 2 

, \ s 4 / ra M 2 m u m d _ m d 2 3 (g 2 GG) m u (qq) m d (qq) 2 

P55l«J - 61 4 4 o 7r 6 +( - 1536tt 6 1536tt 6 15367r eJS 1 6144tt 6 64tt 4 64tt 4 jS 

ml(g 2 GG) _ m u m d (g 2 GG) _ m 2 (g 2 GG) m u (gqaGq) m d (gqaGq) _ (qq) 2 
[ 1024tt 6 1024tt 6 1024tt 6 64tt 4 64tt 4 12tt 2 ,S 

\lmj(qq) 2 7m u m d (qq) 2 17m 2 (gg) 2 m u (g 2 GG) (qq) m d (g 2 GG)(qq) _ (qq) (gqaGq) 
48tt 2 12tt 2 48tt 2 256tt 4 256tt 4 12tt 2 

I \ - ! m " 2 m um d m d 2 3 m u (qq) m d (qq) 2 
Pi2(sj - l 30727r6 + 1536 ^ 6 + m27r6 > s +1 4g ^4 48^ 4 jS 

m u (gqaGq) m d (gqaGq) (qq) 2 5m 2 u (qq) 2 m u m d (qq) 2 5m 2 d (qq) 2 (qq) (gqaGq) 
[ 32tt 4 32tt 4 6tt 2 >S 12tt 2 2tt 2 12tt 2 6tt 2 

( g 2 GG) 2 3m 2 (g 2 GG) 3m 2 (g 2 GG) m M (g 2 GG)(g g ) m d (g 2 GG)(gg) 
Pl3[S > ' 1024tt 6S 1 512^ 6 + 512^ 6 >S 128tt 4 128^ 4 



= 3m 2 (g 2 GG) 3m u m d (g 2 GG) 3m 2 (g 2 GG) _ m u (g 2 GG) (qq) _ m d (g 2 GG) (gg) 
Pl4[S) ( 4096tt 6 2048tt 6 4096tt 6 jS 128tt 4 128tt 4 
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9m 2 (g 2 GG) 9m u m d (g 2 GG) 9m 2 (g 2 GG) 3m u (g 2 GG)(gq) 3m d (g 2 GG)(gg) 
P23{s) { 20487r6 )s+ g4 ^ 4 + g4 ^ 4 , (M) 

(g 2 GG) 2 3m 2 (g 2 GG) 3m 2 (g 2 GG) m M (g 2 GG)(gg) m d (g 2 GG)(gg) 
P24(S) = 1024^ S + ( 512^ 512^ )S + 128^ + 128^ ' ( 16 ) 

m 2 u (g 2 GG) m u m d (g 2 GG) m 2 (g 2 GG) m u (g 2 GG)(gq) m d (g 2 GG)(qg) 
P25[S> 1 4096tt 6 2048tt 6 4096tt 6 jS 384tt 4 384tt 4 ' 1 J 

/ \ / ™« 2 m u m d m d 2 3 ,m u (qq) m d (qq) 2 
P34(S) = ( ^ " 128/ ~ 256^ )S + + "l^ )s (18) 

15m 2 (. 9 2 GG) _ 15m u m d (g 2 GG) _ 15m 2 (g 2 GG) 3m u (gqaGg) 3m d (gqaGg) _ 2(gg) 2 
( 2048tt 6 1024tt 6 2048tt 6 8tt 4 8tt 4 tt 2 )S 

5m 2 u (gg) 2 6m u m d (gg} 2 5m 2 (gg} 2 5m u (g 2 GG) (qq) hm d (g 2 GG)(qq) _ 2(qq) (gqoGq) 

TT 2 TT 2 TT 2 647T 4 64tT 4 TT 2 

_ (g 2 GG) 2 , 3m 2 (g 2 GG) 3m 2 (g 2 GG) m M (g 2 GG)(gg) m rf (g 2 GG)(gg) 

P35lSj - 10247r 6 S +1 512 ^ 6 + 512 ^ 6 J S 128?r 4 128 ^4 ' l 1J J 



III. NUMERICAL ANALYSIS 

To perform the numerical analysis, we use the values for all the condensates from Refs. [25l - l30j : 

(qq) = -(0.240 GcV) 3 , 
(ss) = -(0.8 ± 0.1) x (0.240 GeV) 3 , 
(g 2 s GG) = (0.48 ± 0.14) GeV 4 , 
m u = 5.3 MeV , m d = 9.4 MeV , 

m s (l GeV) = 125 ± 20 MeV , (20) 
(gsqaGq) = -M 2 x (qq) , 
M% = (0.8 ± 0.2) GcV 2 . 

As usual we assume the vacuum saturation for higher dimensional operators such as (0\qqqq\Q) ~ (0|gg|0)(0|^|0). 
There is a minus sign in the definition of the mixed condensate (g s qaGq), which is different with some other QCD 
sum rule calculation. This is just because the definition of coupling constant g s is different f25l [3l|. 

Altogether we took randomly chosen 50 sets of U and 9i. Some of these sets of numbers lead to negative spectral 
densities in the low energy region of interest, which should be, however, positive from their definition. This is due 
to several reasons. One reason is that the convergence of OPE may not be achieved yet for those currents for the 
tctraquark state. Another reason is that some currents may not couple to the physical states properly. Except them, 
there are fifteen sets which lead to positive spectral densities. We show these fifteen sets of ti and 9i in Table HI and 
label them as (01), (02), • • • , (15). They are sorted by the fourth column "Pole Contribution" (PC): 

f s ° e ~ s/M vp(s)ds 

Pole Contribution = ^ „„ 2 . (21) 

f o °° e- s / M lp(s)ds V ' 

The pole contribution (PC) is an important quantity to check the validity of the QCD sum rule analysis. Here, p(s) 
denotes the spectral function. It depends on the ten mixed parameters as well as Mb and sq. We note that tt-7t 
continuum which we shall study later is not included in the pole contribution. By fixing so = 1 GeV 2 , we show the 
PC values in Table U for the fifteen sets. "PC(0.5)", "PC(0.8)" and "PC(1.2)" denote pole contribution by setting 
Mj> = 0.5 GeV 2 , 0.8 GeV 2 and 1.2 GeV 2 , respectively. We find that the pole contribution decreases very rapidly as 
the Borel Mass increases. Since we have discussed the convergence of OPE in the previous section, and found that 
the Dim=10 and Dim=12 terms are much smaller than the Dim=6 and Dim=8 terms, and so it is only the pole 
contribution which gives a upper limitation on the Borel Mass. The Borel window is wider for the former parameter 
sets (1), (2), ■ • ■ , and narrower for the latter ones. It almost disappears for the set (15), whose mass prediction is 
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TABLE I: Values for parameters tt, 0t, the mass range M a , the pole contribution (PC) and the continuum amplitude a(U, 0i). 
The meaning of these quantities are given in the text. There are altogether fifteen sets, which are sorted by the fourth column 
"PC". "PC(0.5)'\ "PC(0.8)" and "PC(1.2)" denote pole contribution by setting M% = 0.5 GeV 2 , 0.8 GeV 2 and 1.2 GeV 2 , 
respectively. 



No 


tl t2 ts ti th 


01 02 9 3 04 05 


M CT (MeV) 


PC(0.5) PC(0.8) PC(1.2) 


a (GeV 4 ) 


(1) 


0.03 0.03 0.73 0.37 0.24 


2.7 3.4 4.7 5.5 3.6 


510 ~ 580 


92% 52% 13% 


1.2 x 10~ Y 


(2) 


0.03 0.92 0.75 0.70 0.03 


5.6 0.80 4.1 2.9 2.5 


510 ~ 590 


90% 46% 11% 


5.5 x KT' 


(3) 


0.25 0.79 0.16 0.95 0.22 


1.8 1.2 6.1 0.44 1.8 


510 ~ 600 


87% 44% 11% 


3.6 x 10" 


(4) 


0.53 0.26 0.93 0.24 0.76 


2.9 0.40 2.0 2.5 3.3 


510 ~ 610 


85% 41% 10% 


1.7 x 10~ b 


(5) 


0.74 0.54 0.74 0.65 0.67 


0.15 3.1 1.4 2.7 6.1 


520 ~ 640 


81% 36% 8% 


1.9 x 10~ b 


(6) 


0.98 0.50 0.12 0.33 0.03 


2.0 4.0 6.3 1.3 1.6 


510 ~ 590 


82% 32% 6% 


5.8 x 10~ s 


(7) 


0.98 0.42 0.84 0.82 0.72 


0.095 1.5 3.7 2.4 3.0 


540 ~ 700 


70% 26% 6% 


4.2 x 10~ b 


(8) 


0.48 0.68 0.58 0.96 0.04 


1.8 2.5 3.0 4.3 3.7 


530 ~ 690 


70% 25% 6% 


1.9 x 10~ b 


(9) 


0.53 1.0 0.99 0.34 0.86 


5.6 4.8 5.3 4.1 0.076 


540 ~ 700 


68% 24% 5% 


4.5 x 10~ b 


(10) 


0.75 0.96 0.32 0.12 0.11 


4.3 2.6 0.93 5.1 2.9 


560 ~ 760 


57% 17% 4% 


9.5 x 10"' 


(11) 


0.31 0.81 0.71 0.10 


4.2 1.8 2.8 5.4 5.1 


570 ~ 780 


55% 17% 4% 


3.2 x 10~ b 


(12) 


0.47 0.40 0.46 0.91 


0.18 1.9 1.9 0.091 0.94 


540 ~ 730 


58% 16% 3% 


2.0 x 10~ Y 


(13) 


0.60 0.26 0.44 0.27 0.24 


3.3 3.6 0.92 5.9 3.7 


620 ~ 850 


43% 13% 3% 


1.7 x 10~ b 


(14) 


0.74 0.73 0.73 0.32 0.28 


1.3 1.3 4.6 3.3 5.6 


620 ~ 850 


42% 12% 3% 


4.3 x 10~ b 


(15) 


0.65 0.55 0.92 0.19 0.96 


4.9 5.2 4.0 5.5 3.3 


730 ~ 930 


25% 7% 2% 


5.4 x 10~ b 



also much different from others. The Borel window should be our working region. However, since the Borel stability 
is always very good when M B > 0.5 GeV 2 , we shall keep the idea of Borel window in mind and work in the region 
0.5 < M| < 2 GeV 2 . On the other side, we shall care more about the threshold value s . 

By using these fifteen sets of numbers, we perform the QCD sum rule analysis. There are two parameters, the Borel 
mass Mb and the threshold value so in the QCD sum rule analysis. We find that the Borel mass stability is usually 
good, but the threshold value stability is not always good. We show the mass range of cr(600), M a , in Table [H where 
the working region is taken to be 0.8 GeV 2 < s < 1.2 GeV 2 and 0.8 GeV 2 < M B < 2 GeV 2 . We find the mass range 
is small when the pole contribution (PC) is large. 

The parameter sets (01)-(06) lead to relatively good threshold value stability. Taking the set (02) as an example, 
we show its spectral density p(s) in Fig [2] as function of s. It is positive definite, and has a small value around s ~ 1.2 
GeV 2 . Therefore, the threshold value dependence is weak around this point, as shown in Fig. [3] for the extracted mass 
as functions of both M B and Sq. We find all the curves are very stable in the region 0.5 GeV 2 < M B < 2 GeV 2 and 
0.6 GeV 2 < s < 1.4 GeV 2 . From the set (02) we can extract the mass of cr(600) around 550 MeV. From other good 
cases, we find that the mass of er(600) is around 550 MeV as well. 




s [GeV 2 ] 

FIG. 2: The spectral density p(s) calculated by the mixed current 77, as a function of s. We show the results of the parameter 
set (02) as an example. 

The parameter sets (07)-(15) lead to the threshold value stability, which is not good. Taking the set (13) as an 
example, we show its spectral density in Fig. 2] as a function of s (left figure), and the extracted mass in Fig. [5] as a 
function of so (upper three curves). The mass increases with so and we cannot extract the mass from this result. 
Many effects contribute to the mass dependence on the threshold value, but for er(600) the ir-ir continuum contribution 
is probably the dominant one. Hence, we add a term p„ n (s) in the spectral function in the phcnomenological side to 
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Borel Mass 2 [GeV' 



FIG. 3: The extracted mass of <j(600) as a tetraquark state calculated by the mixed current rj, as functions of the Borel mass 
Mb and the threshold value so. We show the results of the parameter set (02) as an example. At the left panel, the solid, 
short-dashed and long-dashed curves are obtained by setting so = 0.8, 1 and 1.2 GeV 2 , respectively. At the right panel, the 
solid and dashed curves are obtained by setting M B = 0.5, 1 and 2 GeV 2 , respectively. 




FIG. 4: The spectral density p(s) calculated by the mixed current r), as a function of s. We show the results of the parameter 
set (13) as an example. The left figure shows the full spectral density as given on the left hand side of Eq. (|22|l . while the right 
figure is the one with pT„ T {s) subtracted. 



describe the tt-tt continuum: 

P(s) = fyS(s- M$)+ p„ n (s)+ Pcont . (22) 
where p CO nt is the standard expression of the continuum contribution except the tt-tt continuum. To find an expression 




0.6 0.8 1 1.2 1.4 
s [GeV 2 ] 

FIG. 5: The extracted mass of <r(600) as a tetraquark state calculated by the mixed current rj, as functions of the threshold 
value sq. We choose the parameter set (13) as an example. The solid, short-dashed and long-dashed curves are obtained by 
setting M| = 0.5, 1 and 2 GeV 2 , respectively. The upper three curves are obtained without adding the contribution of the 
tt-tt continuum in the spectral density in the phenomenological side, while the lower three curves are obtained after adding the 
contribution of the tt-tt continuum. 
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for p n „(s), we introduce a coupling 

A ww = (OI^Itt+tt-) . (23) 
The correlation function of the tt-tt continuum is 

IWP 2 ) = * I , , * 5 , ■ -= ( 24 ) 



(27r) 2 (p + q) 2 — m 2 + ie q 2 — m 2 + ie 
and the spectral density of the 7r-7r continuum is just its imaginary part 



p™(s) = Imll^s) = ^lA^I 2 . (25) 

lOTI" 1 V S 

We may calculate A W7r by using the method of current algebra if we know the property of the resonance state. However, 
this is not the topic of this paper. Moreover, in this paper we use a general local tetraquark current to test the full 
space of local tetraquark currents, so we again make some try and error tests, and find that the following function 
leads to a reasonable QCD sum rule result, X n7r ~ s. Hence, we take the spectral density of the tt-tt continuum as 



4 777 

p^{s) = a(t i ,e i ydl-—^. (26) 

We add the continuum contribution p- K - K (s) in the phcnomcnological side and perform the QCD sum rule analysis. 
The values of parameter a(ti,9i) are listed in Table |U After adding the continuum contribution, the threshold value 
stability becomes much better. Still taking the set (13) as an example, we show its spectral density in Fig. |4] as a 
function of s (right figure), and the extracted mass in Fig. [5] as functions of so (lower three curves). We see that now 
the spectral density has a small value around s ~ 1.1 GeV 2 , and the stability of the threshold value is significantly 
improved. 

Hence, we made the same analysis for all the other cases. We found all the cases are good except one, which is 
the case (15), where we are not able to get the desired stability as a function of sq. The mass function has a small 
stability region and increases rapidly with sq. Hence, we consider this case is between the good case and bad case, 
and remove it from the further analysis in this paper. We show several results out of all the good cases in Fig. HI 
which are obtained by using the parameter sets (01), (03), (06), (09), (12) and (14). Wc list the used a(U, 9%) in Table 
1 for all the cases. All the masses behave very nicely as functions of the Borel mass and sq as shown in Fig. [5] In our 
working region 0.8 GeV 2 < sq < 1.2 GeV 2 and 0.8 GeV 2 < M§ < 2 GeV 2 , all the cases lead to a mass within the 
region 495 MeV~ 570 MeV. From this mass range, the mass of a(600) is extracted to be 530 MeV ± 40 MeV. 



IV. THE EFFECT OF FINITE DECAY WIDTH 



After the sq stability has been improved, we notice now that the mass increases systematically with the Borel mass 
as seen in Fig. [5] in all the cases. We therefore try to consider a possible reason of this systematic result. The cr(600) 
meson has a large decay width. We parametrize it by a Gaussian distribution instead of the 5-function for the <r(600). 



P FDW (s) = exp ( - {V ~ S :f X? ) ■ 



The Gaussian width ax is related to the Breit-Wigncr decay width T by ax = T/2.4. We set ax = 200 MeV, and 
Mx = 550 MeV, and calculate the following "mass" : 

M\M B ,s ) = { t r^J ■ (28) 

j e exp ^ ) 

Wc find that the obtained mass M is not just 550 MeV, but increases as Mg increases as shown in Fig. [71 Hence, the 
extracted mass in the QCD sum rule analysis ought to depend on the Borel mass. The amount of the change of the 
extracted mass in the QCD sum rule analysis is similar to the one found in this model calculation. Moreover, we find 
that the finite decay width does not change the final result significantly, which we have also noticed in our previous 
paper 0. 
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FIG. 6: The extracted mass of <j(600) as a tetraquark state calculated by the mixed currents rj, as functions of the threshold 
value so. We choose the parameter sets (01), (03), (06), (09), (12) and (14). The results are shown in sequence. The solid, 
short-dashed and long-dashed curves are obtained by setting M% — 0.5, 1 and 2 GeV 2 , respectively. 
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FIG. 7: The extracted "mass" considering a finite decay width. The solid, short-dashed and long-dashed curves are obtained 
by setting M| = 0.5, 1 and 2 GeV 2 , respectively. 



V. SUMMARY 



In summary, we have studied the light scalar meson er(600) in the QCD sum rule. We have used general local 
tetraquark currents which are linear combinations of five independent local ones. This describes the full space of local 
tetraquark currents which can couple to a (600) either strongly or weakly. We find some cases where the stability of 
the Borel mass and threshold value is both good, while in some cases the threshold value stability is not so good. 
The resonance mass has an increasing trend with so, which indicates a continuum contribution. Hence, we have 
introduced a contribution from the 7r-7r continuum, and obtained a good threshold value stability. The mass of er(600) 
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is extracted to be 530 MeV ± 40McV. Very interesting observation is that the mass increases slightly with the Borcl 
mass. We have made a model calculation by taking the Gaussian width of 200 MeV centered at 550 MeV and try to 
make a sum rule analysis. We see a similar increase trend as seen in the QCD sum rule analysis. 

The continuum contribution exists in the background of the cr(600) meson and it is very important to consider this 
fact in the QCD sum rule analysis for exotic states. We have seen clear tendency of the mass increase with the Borcl 
mass after getting good signal of the threshold dependence. The decay width of cr(600) is related to this increase 
tendency. We are now trying to calculate this by using the three-point correlation function within the QCD sum rule 
approach. The present analysis is very encouraging to apply the QCD sum rule including the continuum states for 
other scalar mesons. Moreover, the continuum contribution should be important in many other resonances such as 
A(1405) etc, which lies in some continuum background. In the future, we will use the QCD sum rule analysis with 
continuum to study various resonances. 
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